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ABSTRACT 

We report on observations of the stellar populations in twelve fields spanning the region between 
the Magellanic Clouds, made with the Mosaic-II camera on the 4-meter telescope at the Cerro-Tololo 
Inter- American Observatory. The two main goals of the observations are to characterize the young 
stellar population (which presumably formed in situ in the Bridge and therefore represents the nearest 
stellar population formed from tidal debris), and to search for an older stellar component (which 
would have been stripped from either Cloud as stars, by the same tidal forces which formed the 
gaseous Bridge). We determine the star- formation history of the young inter-Cloud population, which 
provides a constraint on the timing of the gravitational interaction which formed the Bridge. We do 
not detect an older stellar population belonging to the Bridge in any of our fields, implying that the 
material that was stripped from the Clouds to form the Magellanic Bridge was very nearly a pure gas. 
Subject headings: galaxies: evolution — galaxies: stellar content — galaxies: Magellanic Clouds — 
galaxies: interactions 
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1. INTRODUCTION 

Interactions are known to be an important driver of 
galaxy evolution, but a detailed understanding of their 
influence remains elusive. The Magellanic Clouds are a 
particularly compelling target for investigating the ef- 
fects of minor so-called "harassment" interactions, due 
to their proximity to the Milky Way, their close associa- 
tion over at least the past several Gyr, and their abun- 
dant gas reservoirs, which allow for ongoing star for- 
mation. The strongest evidence that their interaction 
has played an important role in driving the evolution of 
the Clouds lies in the extra-tidal features of the Magel- 
lanic Stream and Magellanic Bridge. Unlike the Mag- 
ellanic Stream, which appears to be a pure-gas feature 
( Guhatha kurta fc R citzcl 1998|), there is a known stellar 
population associated with the Magellanic Bridge, and 
by measuring the ages, chemical abundances, and kine- 
matics of these stars, we can obtain strong constraints on 
the evolution of the dynamical event which formed the 
Bridge, and study in detail how star formation proceeds 
in the wake of such an event. 

The Magel lanic Bridge wa s first reported in H I ob- 
servations bv lHindman et all (|1963l ). and a young stellar 
compo nent (the "inter- C loud population") was discov- 
ered bv llrwin et al.l (11985f). who estimated the age of the 
stars at about 10 8 yr. iDemers fc Battinellil {l998) pro- 
vided the most comprehensive study of the inter-Cloud 
population to date; they observed five fields in the west- 
ern Bridge, and found stars as young as 10-25 Myr in 
both clusters and in a diffuse field component up to 9° 
from the SMC. The young inter-Cloud population prob- 
ably formed in situ in the Bridge, in the wake of the 
Bridge-forming event, making it the nearest example of 
a stellar population whose formation was unambiguously 
triggered by a tidal interaction. Yet surprisingly, no de- 
tailed analysis of the star formation history of these stars 
has been performed to date, and we do not even know 
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the full extent of the population, since no fields along the 
ridgeline of the H I gas have been observed east of the 
midpoint between the Clouds. 

Furthermore, no study to date has specifically searched 
for an older component of the inter-Cloud population, 
which would represent a population of stars that was 
stripped from the Clouds during the Bridge-forming 
event. Tidal forces during an interaction should affect 
both gas and stars, so the inter-Cloud population should 
have an old component, other things being equal. 

The present study will address these open questions 
regarding the inter-Cloud population in the Magellanic 
Bridge. We present our observations and data reduction 
in Section^ In Section I3TTI we trace the eastward extent 
of the young inter-Cloud population, and in Section 13.21 
we search for tidally-stripped stars in our Bridge fields. 
We briefly examine the outer structure of the LMC using 
our four fields nearest that galaxy in Section [3T3l Finally, 
we present the star-formation history of the young inter- 
Cloud population in Section [331 an d summarize the re- 
sults in Section 0J 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. The Observations 

The data were obtained on the nights of January 4 
and 5 2006 (UT), at the Cerro-Tololo Inter- American 
Observatory (CTIO) 4-meter telescope. We used the 
Mosaic-II camera, which images a 36'x36' field onto a 
8kx8k CCD detector array, to obtain short and long ex- 
posures in Washington C, Harris R, and Cousins / filters 
at twelve field positions spanning the inter-Cloud region, 
and at one offset field at a similar Galactic latitude, but 
to the west of the SMC (see Tableland Figure [J). 

The field positions were selected to uniformly sam- 
ple the inter-Cloud region, approximately following the 
ridgelin e of the H I gas which forms the Magellanic 
Bridge (jPutmanl [2000) . In addition, we selected two 
fields to lie off the main ridgeline, but still in regions 
of abundant H I emission (fields mb03 and mbl4). Our 
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fields include the eastern half of the inter-Cloud region, 
where there are very few known star clusters, and where 
iDemers fc Battinellil 's fields mostly lie far from the H I 
ridgeline. 

The exposure times listed in Table [1] were chosen in or- 
der to detect stars as faint as the ancient main-sequence 
turn-off with S/N> 10 in all three filters. At each field 
position and for each of the C, R, and I filters, we 
obtained a pair of long exposures for cosmic-ray rejec- 
tion, and a short exposure to record the photometry of 
brighter stars that are saturated in the long exposures. 
Seeing was stable during both nights; the FWHM varied 
between 0.7" and 1.0". We observed standard star fields 
several times per night for photometric calibration (see 
Table [2]), twilight flats were obtained during evening and 
morning twilight of both nights, and bias frames were 
obtained in the afternoon prior to both nights. 

2.2. Data Reduction 

The data reduction followed the procedure docu- 
mented by the NO A O Deep Wide Field Survey team 
(|Jannuzi et al.|[2003f ). and utilized the mscred package 
in IRAF 1 . Before reducing the data, we obtained an up- 
dated world coordinate system database for the CTIO 
4-meter telescope from the CTIO website, dated from 
May 2004, and we also obtained an updated crosstalk- 
correction parameter file. 

We first processed the bias frames using cedproc and 
zerocombine. We then used cedproc to perform over- 
scan correction and trimming, bias subtraction, bad pixel 
masking, amplifier merging and crosstalk correction on 
all data images. 

We ran obj masks on each twilight-sky exposure to 
mask out any detected stars in the images, and then com- 
bined the exposures using sf latcombine. This yielded 
normalized C, R, and / twilight-sky flats for each observ- 
ing night. We then divided each data frame by the appro- 
priate normalized flat, using cedproc with the sf latcor 
parameter activated. To further flatten the data images, 
we also constructed night-sky flats from the collection of 
data images, after running obj masks on each image to 
mask out detected objects. We then ran cedproc once 
again to apply the night-sky flats. 

We empirically determined the world coordinate sys- 
tem (WCS) for each exposure using msccmatch, which 
uses an astrometric database to identify stars in the im- 
age and determine the WCS from their known coordi- 
nates. We chose to use the Guide Star Catalog, version 
2 for the astrometric database, which we found yielded 
fit residuals that were smaller by about a factor of two 
compared to the USNO-A2 catalog, msccmatch provides 
an interactive interface with which the fit can be im- 
proved by eliminating outlying mismatched points. We 
were able to get the fit residuals in each image below 
0.25 arcsec. 

Next, we used mscstack to combine each pair of long 
exposures into a single, averaged image. We had initially 
activated the cosmic-ray rejection feature of mscstack, 
but we found that if the seeing changed between the two 

1 IRAF is the Image Reduction and Analysis Facility, distributed 
by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foun- 
dation. 



exposures, the stellar flux in the averaged image would 
be significantly clipped by the rejection algorithm. In- 
stead, we stacked the images without rejection, and used 
craverage to detect cosmic-ray hits in the combined im- 
age and replace the affected pixels with the average of the 
surrounding pixel values. We also ran mscstack on the 
short-exposure images, despite the fact that they did not 
have multiple exposures to combine, because mscstack 
also performs a pixel-value replacement for the pixels in 
the bad pixel masks. 

The final step in the reduction pipeline is to run 
mscpixarea, which corrects each image for the variable 
pixel scale across the field. We then use mscsplit to 
separate the eight CCD images per field into separate 
FITS files. We choose to analyze the CCD images sepa- 
rately, in order to properly account for slight differences 
between the CCDs. 

2.3. Instrumental Photometry 

We used the daophot package in IRAF to perform stel- 
lar photometry on our images, using the method of point- 
spread function (PSF) fitting. For each CCD image, we 
ran daof ind to identify peaks, and phot to obtain pre- 
liminary relative photometry of the detected sources. 

Next, we used pstselect to select fifty bright, rela- 
tively isolated sources to form a high-quality sample from 
which an empirical PSF model will be built. Candidate 
PSF stars are rejected if they contain saturated pixels, or 
if their centers are within ten pixels of any pixel flagged 
as bad by our data reduction procedure. 

The PSF fitting is performed in a fully automated 
way, without user intervention. However, we do examine 
residual images after PSF subtraction to ensure that the 
automatic model parameters are correct. All available 
forms for the PSF are explored to find the best-fitting 
model. Next, group and nstar are run to identify and 
photometer stars that are close neighbors of the PSF 
sample. Substar is used to subtract these close neigh- 
bors from the image. We then repeat the PSF fitting, this 
time using the neighbor-subtracted image to produce a 
more accurate model. 

Once the second-pass PSF has been determined, we 
run allstar to perform iterative PSF photometry of all 
detected sources in the image. Then we run daof ind 
on an image in which all known objects have been sub- 
tracted using the PSF model, in order to find fainter stars 
in the image. Then allstar is run on the subtracted im- 
age, using this new list of fainter stars 2 . The two allstar 
photometry lists are combined into a single photometry 
table for the frame, and the IRAF task wesctran is used 
to convert the stars' X,Y pixel coordinates to right ascen- 
sion and declination, using the world coordinate solution 
we determined during the data reduction procedure. 

The above photometry procedure results in some spu- 
rious detections due to two artifacts: bleed columns, and 
scattered-light halos around extremely bright stars. The 
bleed columns are flagged as bad pixels by the data re- 
duction procedure, but this does not prevent daof ind 
from identifying sources along the bleed columns, nor 

2 Alternatively, we could have performed the second-pass allstar 
on the original image using a concatenated list of all detected ob- 
jects. We tested both methods on an image from mb20, the most 
crowded field. The differences in the photometry are consistent 
with the uncertainties estimated by allstar 
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does it prevent allstar from trying to photometer these 
false sources. We therefore remove objects from the pho- 
tometry table that are within 10 pixels of a flagged bad 
pixel. Extremely bright stars have large scattered-light 
halos in the images, resulting in circular concentrations 
of false detections centered on these stars due to the el- 
evated flux levels. To clean the photometry tables of 
sources detected in the wings of extremely bright stars, 
we reject all sources with anomalously high estimated sky 
values. In the absence of these extremely bright stars, 
the sky levels are quite stable, making the sky levels an 
efficient way to identify spurious objects in the wings of 
bright objects. 

The allstar program provides very accurate relative 
photometry, but because the PSF models are uncertain 
at large radii, the allstar magnitudes are normalized to 
a relatively small aperture size of 4 pixels. We there- 
fore need to apply an aperture correction to convert the 
allstar magnitudes to true instrumental magnitudes that 
represent the total flux recorded from each star. To de- 
termine the aperture correction, we select bright, isolated 
objects from our sample and perform concentric-aperture 
photometry using phot with a series of aperture sizes, up 
to 17 pixels. The per-star aperture correction is simply 
the difference between the star's allstar magnitude and 
its large-aperture magnitude: (m Q / s — m a pn). However, 
since individual objects suffer from measurement errors 
and contamination from neighboring objects, we need 
to statistically determine the characteristic aperture cor- 
rection for the entire frame. The distribution of per-star 
aperture-correction values is typically a Gaussian with an 
asymmetric tail to negative values. The Gaussian spread 
is due to measurement uncertainties, and the asymmet- 
ric tail is due to flux contamination, which is never com- 
pletely mitigated by our selection of isolated objects. 

An accurate determination of the aperture correction's 
value and uncertainty requires that we attempt to iso- 
late the underlying Gaussian shape from the asymmet- 
ric skew caused by contaminating flux. To do this, we 
first determine the approximate position of the distribu- 
tion's peak, and then fit a Gaussian function to the points 
to the positive side of this peak value, thereby ignoring 
the negative half that may suffer from flux contamina- 
tion. We adopt the central value of the fitted Gaussian 
as the frame's aperture correction, and its width as the 
uncertainty in the frame's aperture correction. The in- 
strumental magnitude of each star is simply its allstar 
magnitude plus the frame's aperture correction; we also 
add the aperture correction uncertainty in quadrature to 
each star's photometric uncertainty. 

2.4. Standard Star Observations and Photometric 
Calibration 

The final step in our determination of the stellar pho- 
tometry is to place the instrumental magnitudes we have 
measured onto a standard photometric system, using 
standard-star observations. Standard star fields were im- 
aged several times on both of the observing nights. One 
standard field (SA 101) was observed at two separate 
visits on each night in order to measure the effect of at- 
mospheric extinction. 

The standard star observations are presented in Ta- 
ble [H W e selected well-k nown standard fields, first mea- 
sured by iLandoltl (|1973l ), expanded for wide-field CCD 



instruments by IStetsonl (l2000f). and c alibrated for the 
Washington C filter by IGeislerl (|1996l ). The standard 
fields were reduced using the pipeline procedure de- 
scribed in Section [2~2l We then identified sources in each 
standard field with daof ind, and performed concentric- 
aperture photometry on all sources with phot in IRAF's 
daophot package. 

We use the measured aperture photometry of the ob- 
served standard star s, togeth e r with thei r total photom - 
etry as published by IStetsonl (|2000D and IGeislerl (|1996f) . 
to solve the following photometric calibration equation 
for each filter and each CCD in the Mosaic-II array, and 
independently for the two nights of observing: 

AI = m + A + B*(R-I)-C*X 

where M is the published total magnitude, m is the 
observed instrumental aperture magnitude, A is the pho- 
tometric zeropoint, B is the color term, R — I is the star's 
true color, C is the atmospheric extinction term, and X 
is the airmass. Once we have determined A, B and C 
for each CCD and filter, we will be able to convert the 
observed photometry of any star to its total photometry, 
given its observed color and the airmass at which it was 
observed. 

We first determine C, the atmospheric extinction term, 
by examining the photometry from SA 101, the standard 
field that was observed at different airmasses during each 
night. The C parameter is independent of the color term 
and zeropoint, so we simply need to fit a linear regression 
through the observed stars' magnitudes as a function of 
the observed airmass. The slope of the linear regression 
is C, the atmospheric extinction term. Note that we 
need not restrict ourselves to the actual standard stars 
for this step; since we only need the relative photometry 
to determine the extinction term, all of the stars observed 
in field SA 101 can be employed. 

Having determined C, we proceed to simultaneously 
determine the zeropoint and color term. The observed 
magnitudes of all standard stars are first corrected for 
atmospheric extinction (m x = m — C * X); we then 
construct the quantity M — m Xl the difference between 
the the published total magnitude of the star and its 
extinction-corrected instrumental magnitude, and fit a 
linear regression through the M — m x values as a func- 
tion of the published R — I colors. The slope of this 
regression is the color term B, and its zeropoint is A, the 
photometric zeropoint correction. Table [3] presents the 
photometric calibration parameters for each CCD and 
filter, and for each of the two observing nights. 

Note that in Table O the parameters for the Washing- 
ton C filter are the same in all eight CCDs. The reason 
for this is that there are too few standard stars calibrated 
for the C filter to support an independent determination 
for each CCD (see Table [2]), so we were forced to deter- 
mine an average C calibration for the entire mosaic. 

As noted in Table [3J there were too few R and / stan- 
dards present in some of the CCDs to support an in- 
dependent determination of their color terms and zero- 
points. For these cases, we perform a bootstrap esti- 
mate of the parameter values from those published at the 
CTIO website 3 . We determine the mean offset between 

3 http:/ /www. ctio.noao.edu/mosaic/ZcroPoints. html 
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our determined values of the photometric calibration pa- 
rameters, and those published by CTIO, for the CCDs 
that we were able to analyze. We then apply this mean 
offset to the published values of the remaining CCDs, as 
an estimate of what we would have measured if we had 
observed enough standard stars in those CCDs. For the 
uncertainty in these bootstrapped parameters, we simply 
adopt the standard deviation of the mean offset between 
the observed and published values. 

We simultaneously perform a positional match of the 
sources in the C, R and I photometry lists, and apply the 
above photometric calibration to produce catalogs with 
total CRI photometry. We then match objects between 
the CRI catalogs from the short and long exposures of 
each field. For the positional matching, we use a max- 
imum match radius of 0.5"; to be retained in the cata- 
log, an object must be detected in the R band, with a 
matching detection in either C or /. For objects which 
are matched between the short and long exposures, we 
adopt the weighted mean photometry in the final cata- 
log; objects present only in the short or long catalog are 
included as well. The final calibrated composite photom- 
etry catalogs for the twelve observed fields in the Magel- 
lanic Bridge (plus the observed offset field) are presented 
in Table [4] and rendered as pairs of Hess diagrams in 
Figure [21 A Hess diagram is a pixelized color-magnitude 
diagram (CMD) in which each pixel value is proportional 
to the number of stars in the region covered by that pixel. 
The CMD of the offset field is presented separately in 
Figure [H 

2.5. Statistical Subtraction of Foreground/Background 
Contamination 

It is clear from comparing Figures [2] and [3] that the 
stellar populations in many of the fields are dominated 
by foreground Galactic (and background extragalactic) 
contamination. In order to study the underlying inter- 
Cloud populations, we need to first perform a statistical 
subtraction of the contaminant foreground/background 
population. In doing so, we will assume that the popu- 
lation observed in the offset field is representative of the 
contaminant population in each Bridge field (a reason- 
able assumption, given the similar Galactic latitude of 
the offset field). 

We proceed by first determining a scaling factor for 
normalizing the number of objects in the offset field to 
the number of contaminant objects in each Bridge field. 
This is necessary to account for variations in the effective 
area covered by each field (which arise from masking out 
regions contaminated by bad pixels, very bright stars and 
bleed trails). The normalization factor is simply the ratio 
of object counts in the target and offset fields, for a se- 
lected subregion of each CMD that is expected to contain 
only contaminant objects. For the C — R CMD, the nor- 
malization region is defined by the criteria C — R > 2.4, 
R < 22.4 - (C- R), and R > 23.8 -2*(C-R); for the 
R — I CMD the criteria are R - 1 > 1.0 and 18 < I < 20. 
The normalization regions are outlined with dashed lines 
in Figure [3l In the C — R CMD, the objects in the nor- 
malization region constitute 12% of the total number of 
objects, while in the R — I CMD, the fraction is 14%. 
The normalization factor computed for each target field 
ranges between 0.65 and 0.95. We multiply the offset 
field's Hess diagrams by these scaling factors, and sub- 



tract them from the target field's Hess diagrams. The re- 
sultant statistically-cleaned Hess diagrams for each field 
are shown in Figure 2] 

The statistical subtraction was generally successful in 
removing a component from each field's CMD that is 
consistent with the contaminant population in the off- 
set field. However, there are some artifacts present that 
bear explanation. Specifically, the faint end of many of 
the R — I CMDs appear to show an oversubtracted con- 
taminant population. This is simply due to the fact that 
the offset field's R — I CMD has a fainter detection limit 
than that in most of the target fields. 

3. ANALYSIS 

3.1. Extent of the Young Inter-Cloud Population 

Young (age < 1 Gyr) stars provide an unambigu- 
ous tracer of the inter-Cloud population, because no 
foreground or background contaminants are expected to 
share the bright, blue region of the CMD with these stars. 
From p revious work on the young inter-Cloud popula- 
tion bv lDemers &: Battinelli (1998j), we expected to ob- 
serve a population of young stars near the SMC, coin- 
cident with the young cluster population cataloged by 
iBica fc Schmitt] lfl995h . However, we did not know how 
far the young population would extend toward the LMC 
along the H I ridgeline. We isolate stellar populations 
younger than 1 Gyr in the C — R CMD, by selecting those 
stars with R < 20 mag and C — R < mag (see dashed 
lines in Figure [2]) . Stars matching these criteria are ab- 
sent in all of our fields east of mb0 9, which corres ponds 
roughly to the eastern extent of the lBica &: Schmittl clus- 
ters. Interestingly, field mb09 is also near the posi- 
tion along the H I Bridge where the gas surface density 
drops to th e critical thresho ld for star formation of 3- 

4 M Q pc- 2 (IE ennicuttlll989f l. which corresponds to the 

5 x 10 21 c m" 2 seen througho ut the eastern Bridge in Fig- 
ure 4a of iBriins et all (|2005h . West of field mb09, the gas 
density is sustained at a level three times higher, and this 
is where star formation has been active in the Bridge. It 
would seem that the same star formation threshold ob- 
served for disk galaxies holds for this tidal debris envi- 
ronment as well. 

3.2. Searching for Tidally- Stripped Stars in the 
Inter- Cloud Region 

It is perhaps not surprising that the young inter-Cloud 
population is confined to those regions where the gas 
density is relatively high, if we accept the hypothesis 
that these stars formed in situ, following the formation 
of the gaseous Bridge by a recent gravitational encounter 
between the Clouds. However, the tidal forces that pre- 
sumably formed the Bridge should have stripped stars 
and gas with equal efficiency, so we expect to observe 
a population o f such tidally-stripped s tars i n the inter- 
Cloud region. lYoshizawa fc Noguchil (|2003h conducted 
detailed numerical modeling of the stars and gas in the 
SMC, as it orbits both the LMC and Milky Way, in an 
attempt to reproduce the broad physical parameters of 
the Magellanic system. In their best-fitting model, there 
is an abundant population of stars in the inter-Cloud re- 
gion which formed in the SMC, and had been ejected 
into the inter-Cloud region by a tidal interaction with 
the LMC. The tidally-stripped stars should have a simi- 
lar age distribution to the stars in the galaxy from which 
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they were stripped (at least for ages prior to the Bridge- 
forming event when their histories diverged). Since the 
stellar populations in both Magellanic Clouds exhibit a 
prominent red giant branch and a "red clump" horizontal 
branch, these bright features serve as ideal tracers of a 
putative stellar population that had been stripped from 
either of the Clouds during the Bridge- forming event. 

While some of the fields in Figure [4] do show red gi- 
ant branch and red clump features, these older popu- 
lations appear to be confined to the fields nearest the 
SMC (fields mb02 and mb03) or the LMC (fields mbl6- 
mb20). Furthermore, the surface density of these tracer 
populations increases sharply as the galactocentric sepa- 
ration of the field decreases, consistent with populations 
that are bound to the LMC and SMC. In Section [3731 we 
will demonstrate that the red giant populations in fields 
mbl6-mb20 are consistent with a plausible exponential 
disk distribution centered on the LMC. For now, we sim- 
ply conclude that the red giant populations in these six 
fields near the SMC and LMC are very likely composed 
of stars bound to each respective galaxy, and are not 
indicative of a tidally-stripped stellar population in the 
Magellanic Bridge. 

The Hess diagrams of the remaining six fields (mb06- 
mbl4) show no red features that can be associated with 
an old inter-Cloud population. However, the strength 
of this non-detection is limited by the presence of the 
contaminant population. To enhance our sensitivity to 
a potentially sparse old stellar population, we construct 
a composite pair of Hess diagrams from these six "true 
Bridge" fields, and perform a new statistical contaminant 
subtraction on the composite population (see left panels 
of Figure [5]). Even in this composite Hess diagram which 
covers more than two square degrees of the inter-Cloud 
region, there is no detectable trace of an underlying red 
giant branch or red clump feature. 

We can place an upper limit on the surface density 
of red giant branch stars in these six "true Bridge" by 
adding an artificial old stellar population at the distance 
of the Magellanic system (to — M = 18.7 mag, interme- 
diate between the two Clouds) to the composite inter- 
Cloud population. The artificial old stellar population is 
drawn from a theoretical isochrone (|Girardi et al.ll2002h 
with Z = 0.002 and log(age) = 10.0, to which we add 
photometric errors consistent with the data. We modu- 
late the number of artificial stars added until a red giant 
branch is marginally detectable (Figure . We conclude 
from this exercise that there are fewer than 1000 red giant 
branch stars at the distance of the Bridge and brighter 
than R — 23 mag in the observed co mposite popula tion. 
By applying a stellar mass function (|Kroupal 120011 ). we 
can convert the upper limit on the number of observed 
red giants to an upper limit on the total stellar mass 
present in a putative old stellar population. However, 
the conversion factor depends on the assumed age of the 
stars, because the fraction of the total stellar population 
that is brighter than R — 23 mag varies with age. For 
a 10 Gyr population, the upper mass limit is 14800 Mq, 
and for a 2.5 Gyr population, the upper mass limit is 
5300 Mq. Thus, the stellar surface mass density in these 
six "true Bridge" fields is < 0.009 M pc~ 2 ; this is more 
than 400 times smaller than the average surface mass 
density of H I in the Magellanic Bridge (4 Mq pc~ 2 , 
converted from the characteristic column density in the 



Bridge, iBriins et al.ll2005f) . There does not appear to be 
any trace of a tidally-stripped stellar population in the 
Magellanic Bridge, at least in these six fields along the 
H I ridgeline. 

One potential caveat in this analysis is that we have 
assumed that the putative tidally stripped stellar pop- 
ulation would be spatially coincident with the gaseous 
Bridge. This need not be the case; if ram-pressure from 
the Milky Way halo has played a si gnificant role in the 
evolu tion of the Magellanic system ( Ma stropietro et all 
2005), then it is possible that the gaseous Bridge is now 
displaced from the region occupied by tidally-stripped 
stars between the Clouds. We investigate this possibility 
using data from the 2-Micron All-Sky Survey (2MASS, 
iSkrutskie et"all[200l . 

Using the Gator web-based database query service 4 at 
the NASA/IPAC Infrared Science Archive, we obtained 
near-infrared JHK photometry from the 2MASS All-Sky 
Point Source Catalog, in two regions (shown as dashed 
boxes in Figure The first 2MASS region (the "full- 
bridge region") was selected to cover all plausible lo- 
cations where a tidally-stripped inter-Cloud population 
might exist. We selected a range in right ascension be- 
tween 2.5^ and 3.5^, because these limits are bracketed 
by fields mb06 and mbl4, which define the edges of the 
"pure bridge" section of our sample, uncontaminated by 
LMC or SMC stars. We selected a very large range in 
declination, from —77° to —69°, to cover all plausible tra- 
jectories of a putative tidally-stripped stellar population. 
The second 2MASS region (the "SW-LMC region") was 
selected as a comparison field that is known to contain an 
old stellar population at the distance of the Magellanic 
system. This region spans A.2 h to 5 h in right ascension, 
and —75° to —74° in declination. It is coincident with 
our fields mbl8, mbl9 and mb20, in which we have ob- 
served an old stellar population associated with the LMC 
(Section l3.2[) . While the full-bridge region covers a solid 
angle ten times larger than that of the SW-LMC region 
(35 square degrees and 3.2 square degrees, respectively), 
the 2MASS catalog contains about the same number of 
stars in both regions (94000 stars in the full-bridge re- 
gion, and 92000 stars in the SW-LMC region), due to the 
larger stellar surface density of the SW-LMC region. 

The 2MASS J - K CMDs for these two regions are 
shown in Figure [S] In the SW-LMC region, there 
is an abundant population of red objects consisting 
of old stars associated wit h the the LMC. Following 
Nikol aev fc Weinberg! (|2000[ ). we identify the various sub- 
populations of these red objects. The bulk of the pop- 
ulation extends in a narrow diagonal sequence from 
J—K=l mag, If =14 mag to J— Jf =1.25 mag, if =11 mag. 
Along this sequence, there is a sharp drop in the density 
of stars around K=12.3 mag; this is the tip of the red 
giant branch. The stars in this sequence brighter than 
K=12.3 mag are oxygen-rich asymptotic giants, while 
the stars which extend redward of J — if =1.25 mag are 
carbon-rich asymptotic giants. In the full-bridge region's 
J — K CMD, there is a small number of stars whose pho- 
tometry is consistent with these features (notably the ~ 6 
red objects around K = 11 mag which may be Carbon 
stars at the Magellanic distance), but considering the 
very large solid angle covered by the full-bridge region, 

4 http:/ /irsa. ipac.caltech.edu/applications/Gator 
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we do not regard these objects as a significant detection 
of an old inter-Cloud population. We can place an up- 
per limit on the number of red giants at the distance 
of the Magellanic system that can remain undetected 
in the 2MASS CMD, using the same synthetic popula- 
tion analysis described above for our optical CMDs. A 
red giant population containing 150 stars brighter than 
K = 14 mag is easily detectable when added to the 
2MASS CMD, which implies that any old inter-Cloud 
stellar population that may be present has a total stellar 
mass no greater than 2 x 10 6 M@. T his is 1% of the tota l 
H I mass in the Magellanic Bridge ([Brims et al.l [2005) . 
However, the area covered by our full-Bridge region is 
about thr ee times smal ler than the area used to define the 
Bridge by Bri ms et all so the true limit from the 2MASS 
data is closer to 3%. Thus, even accounting for the possi- 
bility that a putative tidally-stripped stellar population 
may be displaced from the gaseous Magellanic Bridge, 
we can still conclude that the Bridge material was more 
than 97% gas when the Bridge was formed. 

3.3. The Outer Disk of the Large Magellanic Cloud 

In fields mbl6-mb20 we observe old stellar populations 
that we conclude are bound members of the LMC, based 
on the sharp increase in their surface de nsity with de- 
creasi ng angular separation from the LMC. lGallart et aT] 
(2004) and others have found that the LMC's stellar 
radial profile follows an exponential disk to projected 
radii beyond 7 kpc, with no sign of a break which might 
indicate the onset of a kinematic halo. Fields mbl6- 
mb20 have projected separations from the LMC of be- 
tween 5 kpc and 8.5 k pc; however, when the orien ta- 
tion of the LMC disk ijvan der Marel fc Cionil 1200 If ) is 
taken into account, the in-disk galactocentric distances 
of these fields are between 6 kpc and 10.5 kpc. We use 
the number of stars in the red clump feature as a proxy 
for the stellar surface density, and plot the surface den- 
sity profile in Figure [7] The solid curve represents the 
best-fit exponential-disk model, with a scale length of 
a = 0.98 kpc, and the dotted and dashed curves are ex- 
ponential disk models fit by previous authors, as noted 
in the figure caption. The surface density profile of the 
red clump stars in fields mbl6-mb20 are generally con- 
sistent with previous measurements of the LMC's outer 
exponential disk, but the fact that the profile is some- 
what steeper in this southwestern quadrant is interest- 
ing. A comprehensive survey of the stellar populations 
in the outer LMC is currently underway; we will there- 
fore postpone further discussion of the LMC's structure 
until this survey is completed, when more definitive con- 
clusions can be made. 

3.4. Characterizing the Purely Tidal Stellar Population 
in the Inter- Cloud Region 

We have determined that the stars in the inter-Cloud 
region appear to be exclusively composed of a stellar pop- 
ulation that formed in situ, in the wake of the Bridge- 
forming event (modulo some contribution from stars still 
bound to the LMC and to the SMC, in the observed 
Bridge fields nearest those galaxies). This isolation of 
a tidally-triggered stellar population provides an impor- 
tant opportunity to examine the nature and evolution of 
star formation processes in tidal debris. We measure the 



age distribution of the inter-Cloud population to deter- 
mine when the star formation occurred, and how long it 
lasted. Since these stars presumably formed in the wake 
of the Bridge- forming event, these measurements provide 
an important constraint on the timing of that event. We 
will also look for spatial structure in the age distribu- 
tion, which may provide insights into how star formation 
proceeds when triggered by a gravitational interaction. 

Previous studies of the inter-Cloud population have es- 
timated the age of the youngest st ars present, through 
simp le isochrone fitting (e.g., iDemers fc Battinellil 
fl998h . Here we will perform a more detailed anal- 
ysis, using the StarFIS H star formation history fitter 
(|Harris &: Zarit skv 2001]). This analysis is motivated by 
the clear presence of composite stellar populations in 
some of our fields, and by our goal to constrain the du- 
ration of star-formation activity in the Bridge. 

StarFISH constructs a library of synthetic CMDs, each 
of which represents a model of what the photometric ob- 
servations would yield, if the observed stellar population 
had a single age and a single metallicity. The model pho- 
tometry is derived from theoretical isoch rones; in this 
case we chose the latest Padua isochrones (|Girardi et al.1 
2002). In order to accurately predict the observed pho- 
tometric distribution in the CMDs, the models include 
a distance modulus, a distribution of extinction values, 
and a detailed model of the photometric errors. The dis- 
tance modulus was simply chosen to be that of the SMC, 
18.9 mag, because the young stellar populations in the 
Bridge are near the SMC on the sky. The distribution of 
extinction values is drawn from regions n ear the eastern 
edge of the MCPS SMC extinction map (jZaritskv et all 
2002). For the photometric errors, we employ an ana- 
lytic model that reproduces the error statistics in the ob- 
served fields. While we usually advocate for an empirical 
model based on artificial stars tests, these tests are only 
strictly necessary when the data images are crowded. In 
the present case, even in our field with the highest stellar 
surface density (mb20), we have detected roughly 109000 
stars in 8192 2 pixels, corresponding to a mean separation 
between objects of almost 14 pixels. 

The StarFISH model library provides synthetic CMDs 
for the range of ages and metallicities thought to be 
present in the observed population; in the present case 
we constructed synthetic CMDs for 16 age bins spanning 
ages 10 Myr to 12 Gyr, spaced uniformly in log(age), 
and for three metallicity bins, Z = 0.001, Z = 0.002, 
and Z = 0.004. The best-fit SFH is found by determin- 
ing the combination of amplitudes modulating these syn- 
thetic CMDs which produces a composite model CMD 
that most closely matches the observed CMD. To take 
advantage of the full CRI photometric data set in de- 
termining the SFH, the fit is actually performed on the 
CMD pair: C - R vs. R and R - I vs. I. 

Because the contaminant population dominates many 
of our observed fields, it is important to account for con- 
taminants in the SFH fit. We could have run StarFISH 
on the statistically-cleaned data set, but we instead chose 
to use the observed data set, and simply include the con- 
taminant offset-field population as an additional ampli- 
tude in the model, in addition to the normal set of syn- 
thetic CMDs. The code will then select a multiplicative 
amplitude factor that optimally accounts for the contam- 
inant population, just as it does for each of the synthetic 
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CMDs. 

The star formation histories of our twelve observed 
Magellanic Bridge fields are presented in Figure [51 and 
the results are consistent with the qualitative analysis of 
the CMDs presented in Sections 13.11 and [3.21 Recent star 
formation has occurred only in the fields west of mbll, 
and old stellar populations are confined to the six fields 
nearest the SMC (mb02 and mb03) and LMC (mbl6- 
mb20). In fields mbll, mbl3 and mbl4 there is no trace 
of a stellar population associated with the Magellanic 
system. 

Our StarFISH analysis shows that star formation in 
the Bridge began around 200-300 Myr ago, and this 
measurement provides an important constraint on the 
timing of the Bridge-forming event. In field mb02, the 
field nearest the SMC, we see a prolonged star formation 
episode spanning ages 80-300 Myr, and in mb03 we see 
a slightly shorter episode spanning ages 100-200 Myr. 
In the other fields in which a young stellar population is 
present (mb06-mb09), the star formation rates are much 
lower, making ages and durations more difficult to deter- 
mine reliably. To boost the signal, we construct a com- 
posite population from these three fields, and determine 
the SFH of the composite population (see Figure [5]) . In 
these more easterly fields, we see evidence for two distinct 
episodes of star formation, 160 Myr and 40 Myr ago. It is 
interesting that while star formation was active through- 
out the western Bridge 100-200 Myr ago, the more recent 
episode 40 Myr ago was apparentl y confined to regions 
furthe r from the SMC. We note that lDemers Sz Battinellil 
(1998) also found that the youngest inter-Cloud popula- 
tions are to be found in fields eastward of the SMC Wing. 

The StarFISH solutions indicate that no significant 
star formation occurred in the Bridge more recently than 
40 Myr ago. This conclusion can be confirmed by direct 
inspection of the CMDs in Figure (2) in no region do we 
see a significant number of main sequence stars brighter 
than R = 15 mag, corresponding to the main-sequence 
turn off position of a 40 Myr isochrone at the distance of 
the Magellanic system. 

The conclusion that star formation in the Bridge 
largely ceased around 40 Myr seems to be at odds 
with a variety of previous research that finds ev- 
idence of__sJell_ai_rjopj^tions much younger than 
this. iDemers fc Battinellil (|1998ft used mam-sequence 
isochrone fitting to conclude that the western Bridge 
contains stars as young as 10-25 Myr. iMeaburnl (|1986l ) 
reported the discovery of DEM 171, a large circular Ha 
filament in the Bridge whic h is likely photoionized by one 
or more massive O stars. iBica fc Schrmttl (|1995| ) found 
that some of the clusters in their Bridge catalog have as- 
sociated emiss ion nebulae, again i mplying the presence of 
massive stars. iMizuno et al.l (|2006f) detected cold molec- 
ular clouds in the Bridge, which demonstrates at least 
the potential for ongoing star formation. 

This apparent contradiction can be partly reconciled 
by understanding that we are not claiming there are ab- 
solutely no stars in the Bridge younger than 40 Myr; 
we find that the star-formation rate dropped off around 
40 Myr ago, and has remained consistent with zero since 
then. We also note that our field selection covers the 
HI ridgeline of the Magellanic Bridge uniformly (see Fig- 
ure []}, except for the segment around RA — 2 h , where 



much of the evidence for more recent star formation 
is to be found. These explanatio ns do not reconcile 
our re sult with the conclusions of IDemers fc Battinellil 
(1998), however. They reported the widespread presence 
of stars aged 10-25 Myr in a number of fields east of 
RA = 2 . This is based on an analysis of their Fig- 
ure 7 (lower panel), in which theoretical isochrones are 
overplotted on a composite CMD from four of their ob- 
served fields. While the observed main sequence does 
appear to follow the shape of the 10 Myr isochrone, it 
is clearly truncated around My = — 3 mag, whereas a 
10 Myr population should have a main sequence that ex- 
tends up to My = —5 mag. A main sequence turn-off 
at My = — 3 mag is consistent with the 40 Myr age that 
we have found for the youngest bulk population in the 
Bridge. 

4. SUMMARY 

We have observed stellar populations in twelve fields 
uniformly spanning the region between the Magellanic 
Clouds. Our fields were selected to follow the ridgeline of 
the H I gas that forms the Magellanic Bridge, in order to 
look for stars that formed in situ in the Bridge from gas 
that had already been removed from one of the Clouds, 
and also for stars that were stripped from either of the 
Clouds by the same tidal forces that presumably stripped 
the gas. 

We observed the previously known young stellar popu- 
lation in the western half o f the inter-Cloud region, mos t 
recently characterized by IDemers fc Battinellil <|1998h . 
and extend on previous analyses in two key ways. First, 
we determine that the eastward extent of these stars 
is truncated around a = 3 , corresponding also to 
the eastward extent of the star clusters cataloged by 
IBica fc Schmittl (|1995l ). and to the point at which the H I 
surface density falls below t he critical thresho ld for star 
formation as determined bv iKennicuttl (l989). Second, 
we use the StarFISH program to determine quantitative 
star formation histories of the young inter-Cloud pop- 
ulation, finding that star formation in the Bridge com- 
menced about 200-300 Myr ago, and continued over an 
extended interval, until about 40 Myr ago. 

We found no evidence for a population of tidally- 
stripped stars in the inter-Cloud region, and our non- 
detection allows us to conclude that the material stripped 
from the Clouds into the Bridge was very nearly a pure 
gas, with an upper limit on the mass fraction in stars 
of less than 10~ 4 if coincidence with the gaseous Bridge 
is assumed, and 0.03 otherwise. This can potentially be 
understood if the pre-collision SMC had an extended en- 
velope of gas, surrounding a more tightly bound stellar 
component. In this weak tidal interaction might 

unbind the gas envelope while leaving the stellar compo- 
nent undisturbed. It is known that some dwarf galaxies 
have H I gas extending beyond 2- 3 times the radii oc- 
cupie d by their stellar populations (Sal peter fe Hoffman! 
1996), so perhaps this scenario is plausible. In fact, 
the recent numerical simulation of the tidal history of 
the SMC by citeyn03 included such an extended gas 
envelope, in order to produce a pure-gas Magellanic 
Stream. Nevertheless, the inter-Cloud region in their 
best model contains an abundant populati on of tidally- 
stripp ed stars. In addition, Figure 4a of iBriins et al.l 
(2005) shows that the H I gas in the Magellanic Bridge 



appears to be contiguous with the higher-density gas 
in the central regions of the SMC, which are currently 
abundantly populated with stars. If the gaseous Bridge 
formed via the tidal extraction of this high-density gas 
from the central regions of the SMC, then the question 
remains: where are the stars in the Bridge that should 
have felt these same tidal forces? 

Future observations may be able to address this ques- 
tion. While some kinematic measureme nts exist for a 
hand f ul of stars in the int er-Cloud region aurice et al.1 
Il987t iKunkel et al.|[l997l ). radial- velocity kinematics of 
a truly representative sample of the young inter-Cloud 
population would help us to better understand the dy- 
namical evolution of the Magellanic Bridge. A much 
deeper understanding of the SMC's complex three- 
dimensional structure and kinematics would certainly 
help as well. These measurements (along with our cur- 
rent understanding of the orbital motions of the Clouds 
and the Milky Way) could then be used to motivate new 
detailed numerical simulations specifically targeting the 
formation of the Magellanic Bridge as a pure-gas feature. 



We may then better understand how tidal features are 
formed during minor harassment interactions, and what 
role such interactions play in driving the evolution of the 
participant galaxies. 
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TABLE 2 
Standard Field Observations 
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5 








R 


10 


2006-01-04 00:55 


1.38 




89 








I 


10 


2006-01-04 00:57 


1.39 




132 


SA 92 


00" 55 m 


00° 40 ; 


(j 


20 


^UUU U-L U1 UO. J.O 


3.01 




6 








R 


10 


2006-01-04 03:15 


3.08 




98 








I 


10 


2006-01-04 03-1 7 


3.14 




147 


SA 101 


Ogh 57m 


—00° 20' 


c 


20 


2006-01-04 05:43 


1.34 




6 








R 


10 


2006-01-04 05:45 


1.33 




51 








I 


10 


2006-01-04 05:47 


1.33 




79 


SA 101 


09 h 57 m 


—00° 20' 


Q 


20 


2006-01-04 07-47 


1.15 




6 








R 


10 


2006-01-04 07:50 


1.15 




49 








I 


10 


2006-01-04 07:51 


1.15 




82 


SA 98 


06' 1 52 m 


—00° 24' 


Q 


20 


2006-01-04 07-^5 

iiUUU Ul KJ 1 .<JtJ 


1.73 




15 








R 


10 


2006-01-04 07:57 


1.75 




662 








I 


10 


2006-01-04 07:59 


1.77 




888 


NCJP 99QS 


06^ 49 m 


—36° 00' 


(j 


20 


9f)0fi 01 04. flK-fl9 


1.36 












R 


10 


2006-01-04 08:04 


1.37 




515 








I 


10 


2006-01-04 08:06 


1.38 




42 


SA 92 


00" 55 m 


00° 40' 


C 


20 


2006-01-05 01:01 


1.42 




7 








R 


10 


2006-01-05 01:03 


1.43 




97 








I 


10 


2006-01-05 01:05 


1.44 




144 


SA 101 


09 h 57m 


-00° 20' 


C 


20 


2006-01-05 04:11 


1.91 




6 








R 


10 


2006-01-05 04:13 


1.88 




55 








I 


10 


2006-01-05 04:15 


1.86 




80 


SA 101 


09 h 57 m 


-00° 20' 


C 


20 


2006-01-05 06:23 


1.22 




7 








R 


10 


2006-01-05 06:25 


1.22 




52 








I 


10 


2006-01-05 06:27 


1.22 




77 
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TABLE 3 

Photometric Calibration Parameters 



Filter CCD Extincton Ext. Unc. ZP ZP Unc. Color Term Col. Term Unc. 



Night 1 


6 


l-8 a 


0.275 


0.0015 


0.074 


0200 


-0.100 


0.0100 


R 


l b 


0.072 


0.0010 


-0.690 


0200 


-0.032 


0.0130 


R 


2 


0.072 


0.0010 


-0.706 


0.0045 


-0.024 


0.0077 


R 


3 


0.072 


0.0010 


-0.684 


0.0023 


-0.042 


0.0039 


R 


4 


0.072 


0.0010 


-0.670 


0.0045 


-0.049 


0.0083 


R 


5 b 


0.072 


0.0010 


-0.682 


0.0200 


-0.035 


0.0130 


R 


6 


0.072 


0.0010 


-0.680 


0.0033 


-0.044 


0.0050 


R 


7 


0.072 


0.0010 


-0.666 


0.0027 


-0.052 


0.0042 


R 


8 


0.072 


0.0010 


-0.703 


0.0027 


-0.030 


0.0050 


1 


1 


0.040 


0.0009 


-0.029 


0.0185 


-0.003 


0.0140 


I 


2 


0.040 


0.0009 


0.009 


0.0022 


0.007 


0.0016 


1 


3 


0.040 


0.0009 


0.008 


0.0020 


-0.001 


0.0016 


I 


4 


0.040 


0.0009 


0.014 


0.0038 


0.005 


0.0041 


I 


5 b 


0.040 


0.0009 


-0.006 


0.0230 


0.003 


0.0090 


I 


6 


0.040 


0.0009 


-0.011 


0.0020 


0.012 


0.0015 


1 


7 


0.040 


0.0009 


0.010 


0.0015 


0.003 


0.0013 


I 


8 


0.040 


0.0009 


-0.003 


0.0031 


0.010 


0.0022 


Night 2 


C 


l-8 a 


0.299 


0.0012 


0.026 


0.0200 


-0.087 


0.0100 


R 


l b 


0.090 


0.0010 


-0.731 


0.0200 


-0.021 


0.0130 


R 


2 


0.090 


0.0010 


-0.732 


0.0078 


-0.022 


0.0140 


R 


3 


0.090 


0.0010 


-0.707 


0.0030 


-0.032 


0.0050 


R 


4 


0.090 


0.0010 


-0.727 


0.0088 


-0.047 


0.0176 


R 


5 b 


0.090 


0.0010 


-0.723 


0.0200 


-0.024 


0.0130 


R 


6 


0.090 


0.0010 


-0.710 


0.0055 


-0.030 


0.0082 


R 


7 


0.090 


0.0010 


-0.727 


0.0067 


-0.022 


0.0094 


R 


8 


0.090 


0.0010 


-0.748 


0.0043 


-0.020 


0.0079 


I 


l b 


0.061 


0.0007 


-0.038 


0.0220 


0.007 


0.0130 


I 


2 


0.061 


0.0007 


-0.011 


0.0025 


0.002 


0.0017 


I 


3 


0.061 


0.0007 


-0.030 


0.0047 


0.004 


0.0038 


I 


4 


0.061 


0.0007 


-0.003 


0.0111 


-0.019 


0.0109 


I 


5 b 


0.061 


0.0007 


-0.033 


0.0200 


0.001 


0.0130 


I 


6 


0.061 


0.0007 


-0.029 


0.0026 


0.018 


0.0021 


1 


7 


0.061 


0.0007 


-0.021 


0.0031 


-0.002 


0.0021 


I 


8 


0.061 


0.0007 


-0.029 


0.0028 


0.011 


0.0020 



a There were to few C standards to support an independent transformation solution for each CCD. b There were too few standards in this CCD to 
support an independent transformation solution. The zcropoint and color term values for this CCD are bootstrapped from the values published at 
the CTIO website. 



TABLE 4 

Photometry of Stars in Magellanic Bridge Fields 



Object ID 


Right ascension 


Declination 


C 




R 




I 


°7 








[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Region mb02 


1 


m n 43 m w 7 s 


-74° 37' 51" 


19.187 


0.037 


19.416 


0.028 


19.385 


0.039 


2 


01 h 43 m 10.7 s 


-74° 38' 52" 






21.070 


0.087 


20.796 


0.115 


3 


01 h 43 m W 9 s 


-74° 37' 36" 


18.075 


0.037 


18.486 


0.026 


18.482 


0.032 


4 


Ql h 43m n o s 


-74° 44' 39" 


21.457 


0.116 


20.501 


0.052 






5 


0l h 43 m 11.2 s 


-74° 33' 20" 






20.255 


0.045 


19.720 


0.059 


6 


Q1 h 43m U 2 s 


-74° 36' 26" 


20.416 


0.060 


18.901 


0.023 


18.379 


0.033 


7 


Ql h 43 m n 2 s 


-74° 37' 22" 


20.728 


0.065 


19.807 


0.041 


19.424 


0.038 


8 


Ql h 43 m n 5 s 


-74° 36' 02" 






19.034 


0.030 


18.538 


0.033 


9 


Ql h 43 m n 8 s 


-74° 48' 51" 


21.007 


0.071 


20.693 


0.078 






10 


Ql h 43 m U gs 


-74° 45' 31" 


17.561 


0.035 


15.270 


0.032 







Note. — 



The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
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5 4 3 2 

RA (hours) 



Fig. 1. — A schematic view of the Magellanic Bridge regio n. The positions of the tw elve observed CTIO-4m fields are shown as 
numbered boxes . The inter-Cloud fields previously observed by Demcrs & Battinclli (1998) arc shown as gray boxes. Star clusters from 
Bica & Sc hmltg llT995t) are shown as points. The gray line is an approximate trace of the HI ridgeline through the inter-Cloud region, from 
Figure 1 in Putman (2000). The h ashed areas re present the portions of the regions covered by the Magellanic Clouds Photometric Survey 
(LMC at left. lZaritskv et al.ll200l SMC at right. IZaritskv et al.ll2002h . 
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Fig. 2. — Color-magnitude Hess diagrams of each of the twelve fields observed in the Magellanic Bridge. Each pair of panels consists 
of the C — R,R CMD (left) and the Ft — 1,1 CMD (right). The fields are numbered in order of increasing right ascension (see Figure [TJ. 
Dashed lines indicate the region of the C — R CMDs where main sequence stars younger than 1 Gyr are expected. The pixel values indicate 
the number of stars present in each pixel, ranging from zero (white) to 20 (black), as indicated by the scale bar on the right. 
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Fig. 3. — The CMD pair for the offset field, which has a similar Galactic latitude to the Bridge fields, but lies to the west of the SMC. 
It therefore should represent a pure foreground Galactic and background extragalactic population. As in Figure [2] the C — R,R CMD is 
shown at left, and the R — /,/ CMD is shown at right. The normalization regions, which are expected to contain only contaminant objects, 
are outlined with dashed lines. 
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201 201 201 201 
C-R R-l C-R R-l C-R R-l C-R R-l 



Fig. 4. — Color-magnitude Hess diagrams of the twelve observed Bridge fields, after statistically subtracting a contaminant component 
derived from the offset field population (see Figure[3]l. Darker pixels indicate a local excess of stars compared to the contaminant population, 
whereas lighter pixels indicate a local deficit of stars. In a perfect subtraction, the pixels lighter than the medium gray of zero difference 
would be consistent with the noise of the subtraction. Many of the fields show large deficits of stars at the faint end of the R — I CMD; we 
believe this is due to a mismatch in the faint limit between these target fields and the offset field. The pixel values indicate the difference 
in star counts (Nfi e id — N ff se t) for each pixel, ranging from -5 (white) to +10 (black), as indicated by the scale bar on the right. 
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observed + 500 RGB + 1000 RGB + 2000 RGB +4000 RGB 
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02 02 02 02 02 
C-R C-R C-R C-R C-R 



Fig. 5. — The top left panel shows the C — R Hess diagram for a composite population formed by combining the "true Bridge" fields 
(mb06-mbl4). We see no trace of a red giant branch feature which would indicate an older stellar population in these fields. We determine 
the upper limit on the number of red giants which could remain undetected in this diagram by adding a synthetic red giant branch composed 
of increasing numbers of stars. From left to right, starting with the second panel we show the same Hess diagram after having added 500, 
1000, 2000, and 4000 synthetic red giant stars. The bottom row of panels show the same Hess diagrams, after having statistically subtracted 
a contaminant population derived from the offset field. 
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Fig. 6. — Near-infrared J — K color-magnitude diagrams for two regions between the Magellanic Clouds. At left is the "full-bridge 
region" , which spans 35 square degrees and is designed to cover the full region between the Clouds that is not contaminated by stars bound 
to either galaxy. At right is the "SW-LMC region", which covers 3.2 square degrees and is coincident with our fields mbl8, mbl9, and 
mb20, in which we have observed an old stellar population associated with the LMC. The dotted lines in each panel indicate the region 
where the bright red giant and asymptotic giant branches are expected at the distance of the Magellanic system. While an abundant 
population of these stars is present in the SW-LMC region, the much larger full-bridge region contains hardly any such stars. 
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Fig. 7. — The number of red clump stars present in fields mbl6-mb20 vs. the physical distance of each field from the center of the 
LMC. In determining these distances, we have assumed that the red clump stars are in the disk of the LMC, whose geometry is known. 
The best-fit expone ntial-disk mo del for these data i s shown as the soli d line, and it is similar to previous determinations of the LMC's 
exponential disk by (Hardy 1978, ; dashed line) and (Gallart ct al. 2004, ; dotted line). 
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Fig. 8. — StarFISH SFH solutions for the twelve observed fields in the Magellanic Bridge. Each panel is divided into three subpanels, 
showing the star formation rate as a function of time for three metallicitics (top to bottom: Z=0.004, Z=0.002, and Z=0.001). The best-fit 
history is plotted with a heavy line in each sub-panel, and the shaded gray regions indicate the estimated uncertainty of the solution. Fields 
arc presented in order of increasing right ascension, starting with mb02 (the field nearest the SMC) in the lower left, up to mb20 (the field 
nearest the LMC) in the upper right. 
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Fig. 9. — The SFH solution for a composite population formed from fields mb06-mb09. The solutions for these fields presented in 
Figure[8]are inconclusive, due to the small number of stars in these fields. By forming a composite population, we can better constrain the 
average history of these three regions. 



